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ABSTRACT 

Bathymetric charting and coastal mapping using MSS data 
are two parts that constitute this thesis work. The Landsat 
products such as Computer Compatible Tape (CCT), positive films 
paper prints were obtained from the National Remote Sensing 
T^ency (NASA), Hyderabad, The area chosen for the study is 
Andhra Pradesh coastal region near Mach li pat num. The area is 
covered by one Landsat - 4 scene (Number : 142 : 049), The 
toposheets of this area having the contours of water depth of 
9 m and 18 m, were obtained from Survey of India, Eahradun, 

A relation is developed between reflectance value (Band 
and water depth. This relation has been used in preducting the 
water depth from reflectance value for other areas. Invest igat 

I 

have shown a strong correlation between Band 1 reflectance valu 
(MSS- Landsat 4) and shallow water depths of the sea. The press 
technique enables us to get an up-to-date hydrographic inforraat 
which has been of utmost importance to navigators and coastal 
engineers. The bathymetric charting by Remote Sensing is not 
only time saving but also time effective. 
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Ca^^PTER 1 
INTRODUCTION 


1 .1 GENERAL 

The last few years have seen tremendous growth in the 
field of remote sensing. October 1957 marked the beginning 
of satellite remote sensing with tha launch of SPUTNIK by 
U.S.S.R, For the first time in the history, a rnan-made object; 
v;as circling the earth^, covering different regions and countrii 
of the world from its vantage point in space.. This led to rap: 
development of satellite borne sensor for reconnaissance purpos 
Tlius remote sensing has advanced from visual analysis of aeria? 
photoaraphs to automatic computer analysis of digital multibaj 
imageries. Today, one has nearly lost count of satellites thai 
spy the earth every moment for military and civilian purposes,; 

I 

The civilian remote sensing program was essentially an 
U.S.A, and U.S.S.R. effort in the seventies, U.S.A. made the ; 
data available from its satellites to many countries around th« 
world at extremely reasonable costs. As a consequence the use I 
of such data has ’oecome common around the world, With the 
development of high resolution satellites like FRENCH SPOT j 
system, remote sensing techniques are being reliably applied ir 
large number of fields. A number of studies like crop pattern^ 
determination, crop estimation, landhse, surface water distribij 
river course monitoring, forestry planning, geolo-gical mapping, 
are being taken up. ; 


Here in India, we have various facilities like data 


products laboratory at SPACE APPLICATIONS CENTRE (SAC), full- 
fledged acquisition and processing facility at NRSa, HYD'.LlABAr 
for undertaking various projects for diverse applications. The 
launching of ir-miAN REMOTE SENSING SATELLITE (IIRS-I), which is; 
scheduled by the end of 1937 will give the indigenous efforts a 
further boost in the field of remote sensing. 


1,2 MULTI spectral SCANNER (MSS) 

The MSS oscillating mirror scans the information of the i 
object in cross track direction and it travels along the track : 
as s'nown in Fig, 1,1, It covers an area of 135 kni x 135 km in 
one scene. The information is obtained pixel by pixel and : 
scanned line by line. The total area of 1S5 km x 135 km covere; 
in one scene is divided into 2340 scan lines and each scan line 
is divided into 3240 pixels as shown in Fig, 1,2, but there may 
be small variation in the above numbers. For the data obtained 
for the present investigation the number of scan lines ^^?ere 
2400 and number of pixels as 3233. Pixel is the minimum unit 
that can be distinguished from Landsat data. For MSS, the pixe 
has the dimensions of 57 m x 79 m. 

The MSS sensor scans the data in four different bands of; 
the optical and near IR spectrum simultaneously. The bands ari 
numbered as 4, 5, 6 and 7 in Lands at s 1, 2 and 3 and as 1, 2, 3; 
and 4 in Landsat s 4 and 5, Each one of the band is useful for 




^ultispcctrol scQnncr 
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one particular field of study. The various uses of each 
band are given in Table 1.1. 

1.3 THEMATIC MAPPER (IM) 

It is a special scanner which has better resolution 
than MSS and operates in 7 channels. Various bands are 
numbered as 1 to 7 and each one has application in a specific 
field of study. All the channels have a spectral resolution 
of 30 m against the 79 m in MSS, except for thermal band i.e, 
band 7. The various bands and their uses are given in Table 1.1 

1 . 4 APPLICATION OF REfiOTE SENSING 

Remote Sensing has proved its worth in various 
applications, Indian scientists at various well equipped 
research organisations and frontline academic institutes have 
applied the satellite products fear diverse applications like 
geological, agricultural, landuse with diverse techniques 
like computer aided, visual and image processing. 

Snow melt run-off studies (Ramamoorthi, 1983), ground 
water table prediction (Rampal, 1934), soil mapping (Karale, 
1983), forest survey and management (Madhavan Unni, 1983) 
land evaluation and classification for agriculture (Murthy, ! 
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Venkataratnam and Saxena, 1983) are few to quote anong 
Indian works, 

1*5 CRITERIA FOR TOE SELECTION OF BAND 

It is very important to select a suitable band for 
a given field of study. The criteria for the selection 
are as follows. 

The imageries which are available in Band 1, 2, 3 
and 4 differ in appearance which depend upon the varying 
reflectance of the different features of the earth's surface,! 
Maximum transmission of light in clear sea water is in the 
blue spectral region. Water acts as an optical filter, 
progressively absorbing radiant energy at longer wave lengths; 
until almost complete absorption occurs in the near infrared | 
region. It is shown in Pig, 1,3, 

Shallow clear water or water containing suspended 
materials returns considerable amount of light in less 
attenuated Band 1 which provides valuable information regardli 
the <tepth car nature of foreshcare and bottcana, concentration of 
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suspended material, location of submerged objects. 

The attenuation coefficient is higher in band 2, 3 
and 4 and very little light is directed towards satellite by 
water bodies ii7hiGh appears very dark as compared to land surf a 
Tlierefore, long wavelength bands particularly band 4 is very 
useful in distinguishing land from water and are helpful in 
delineation of coast line, islands, sandspits etc. 

Hence for the present work, band 1 was selected for 
determining water depth and band 4 X'^as selected, for mapping tlj 
land and water. 

1.6 NSSD OF RS^IOTS SSNSIWG FOR HYDROGRAPHIC SURVEYING 

A study of the hydrographic surveying reveals that 
except for the extremely small area (5%) surveyed, the 
remaining areas have been inadequately surveyed, which I 

constitute enormous surveying task. The 1982 Expert Group 
Report on Hydrographic Perspective Plan has assessed that it 
would take a span of 10 years by employing 16 surveying vessel 
and over 4500 personnels to complete navigational surveys anc 
present strength of survey vessels and manpower is deficient. | 
Apart from this, the coastal areas and mouths of rivers on eas 
coast of India like Hugli, Mahanadi, Godavari, Krishna are ! 
changing frequently due to floods and storms. The position, | 
extent and depths of many offshore islands, banks, shoals etci 
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x-/hich are not connected with main land by triangulation are 
required to be determined. Some ships require navigable area 
which are more than 22 m deep. An up-to-date hydrographic 
information has always been of utmost importance to navigators 
and coastal engineers. 

The hydrographers have traditionally relied on sounding ; 
machine^ echo sounding and aerial photography to maintain and 
update ta ./ix: nautical charts. Though these inventory technique 
are quite effective and recent but are slow, cumbersome and 
costly. At the same time catastrophic events like floods, store 
can change the hydrography quickly in delta regions, sandy shori 
etc. Therefore, the hydrographer requires a method which can 
give him quick results on bathymetry, coastline, shoals etc. i 
which he is unable to get by conventional method, without losinc 
much t ime , 

In the continued attempt to introduce advance technology! 
fitting to the navigator's requirement attempts have been made; 
to develop the new nethod of collecting offshore hydrography ; 
through remote sensing techniques. This method has greatly redij 
the cost and improved the quality of work. This has great i 

potential to collect large amount of hydrographic data very ; 

rapidly over shallow coastal water. I 



In this present project, an attempt has been made to 
correlate the reflectance of water bodies obtained by sensor 
and water depth for shallow water. We have prepared the map 
of coastal area distinguishing land, clear water and water 
having suspended materials. Our area of study in Andhra 
Pradesh coastal area between Madras and Kakinada. Details 
of the study are explained in the following chapters. 

1.7 OBJECTIVE OF THE STUDY 

The objective of the study is 

(1) to correlate the reflectance (MSS Band 1) with shallow 
water depth 

(2) to map the coastal area distinguishing land, water 
having suspended materials and clear water. 

To achieve the first objective of the work, it is 
necessary to get water depth contours of that area, to calcui 
the latitudes and longitudes of the points on the known water 
depth contours, to convert these coordinates to line numbers 
and pixel nunabers and to read the corresponding reflectance 
values from CCT obtained from NRSA in MSS Band 1. The 
larocedure for this work is given in chapter 3, 
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For the second objective an area, which covers the 
land, water containing suspended materials and clear water, 
was chosen on imagery and reflectance values of all the pixe] 
falling over that area were read from CCT, Maps were prepare 
by the program which is given in Appendix 6, 
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CHAPTER 2 
SCOPE OP THE WORK 

2*1 EARLIER STUDIES 

Not much literature Is available on the application 
of Remote Sensing in the field of measurement of water depth. 
Bathymetry is very old field for navigators and they are usinc 
so many traditional methods of water depth measurement like 
sounding machine, echosounding, aerial photography etc. Thouc 
these techniques are quite effective but are slow, 

( o \ 

The work of Ross ' and Polcyn and Rollin have 
shown that there is a correlation between MSS band 1 imagery a 
shallow water depth. Correlation was found between MSS - 1 
Imagery and depth in number of areas corresponding to water 
depth of less than 2 m, 5 to lO m. But the depth predicted 
from imagery was approximate. Under favourable condition 
Landsat imagery can detect sub-surface detail as small as 200 
in diameter. Water depth over a large area of world coastline 
and shallow sea can be assessed by Landsat imagery and Image ! 
processing technique which would be valuable aid to chart 
revision by directing the hydrographer to area requiring closd 
examination* | 



The wave refraction method was applied by Polycn in 
1969, The underlying principle in this case is that when a 
wave moves through a liquid the characteristics of the medium 
which may influence its motion are 

(a) The depth and other boundary conditions, 

(b) Gravitation, since the changed profile or contour 
of the surface involves work against gravity. 

(c) Surface tension, because the pressure under a curved 
surface is different from that beneath a flat surface, 

(d) Viscosity, which is the dissipative energy agent. 

The velocity of propagation for a surface wave in a 
liquid of density is given by 

--2 _2^S g"?v 

V = -7^ “lir” ^ 2 . 1 ) 

where S = surface tension 

g = acceleration due to gravity 
Ti <a wave length 

For sufficiently large value of i,e,, for long waves 
first terra is negligible compared to second term of the above 
expression, Tnus the velocity of 'Gravity Waves' is given by 
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V 


2 


* a2i 
2 


These gravity waves are observed over deep part of the water 
body, We know that 


V = f7\ 

where, £ = frequency of wave. 

Since time period of a wave Is the inverse of the frequency 
l,e, time period (p) = — ~ 

Finally we get the relation between velocity and time period as 

V » “ 5,09 p 

Jr 

2 

where g = 32 ft ,/sec 
V being measured in ft ,/sec. 

The velocity in the shallow water is given by the 
expression 

v^ «= gd (2,2) 

where v = velocity of the wave 

g » gracltatlonal acceleration 
d = depth of water. 



As the wave moves from deeper portion to the shallow 
part, the time period of the wave 'p' remains constant. This 
property is utilized to determine the depth. If v^^ be velocity 
for deep water, then we have 



JIT 


(foA©)' 




where = frequency of wave in deep water 

9 ^ 0 ® wave length in deep water 

p = time period of wave = 

Therefore, p^ ® (2,3) 

g 


For shallow water 


2 

Therefore , p 



Now from equations 2,3 and 2,4, we get 


2l 2nd 


and also 



2l 



(2,4) 


(2.5) 



The ratio 


is thus functionally related to the ratio 


A measurement o£ & 7\ will then allow 'd' to be determined. 

This can be done by successive picture at diort Intervals to 

obtain the values of & v. Measurement of can be made 

o o 

directly from an aerial photograjAi. This method can, therefore, 
be applied easily where seashore or the bed of the water body 
slopes gently. 

(3) 

Polcyn & Lezenga demonstrated the use of satellite 
data for mapping of shallow water features for the purpose of 
upgrading the world navigation charts. A mathematical model for 
depth measurement using the ratio of reflectance in Band 1 and 
Band 2 had been successfully developed. Satellite data also 
provide geographical evidence for verifying existence or non- 
existence of doubtful shoal features appearing in the world charts; 
and considered to be hazardous to shipping. 

The technique employs the relation shown by equation 

1 “ 1 ^1 «1 
2 „ « K 2 ) 

Where Z » water depth, 

Kj,K 2 *= attenuation coefficients of water in two different wave 
lengths. 

“ reflectances for bottom material in two different bands] 
“1 “2 * St ants of instrument which are known. 


” incoming solar radiations, 

“ analog signals observed in the multi spectral scanning, 

9 = observation angle. 

0 = solar zenith angle . 

Efere and V 2 are obtained from Band 1 and Band 2 of CCT, 
The incoming solar radiation H 2 ) is available from standard 

references, ^ and are reflectances for bottom materials, 
which are known. 

Photogramme trie bathpaetry is one of the traditional 
techniques of remote bathymetry. Colour aerial photography 
with its remarkable clear water penetration characteristics and 
dramatic presentation of submerged detail is basic tool in this 
method for mapping seabed in water of moderate depth. Aerial 
photogrammetry of bathymetry requires several minor departures 
from the method normally used where bundle of light rays passes 
only through atmosphere. The effect of refraction at the water 
air interface must be taken into account in the aerotriangulation 
for imaged point. The solution of this problem requires a mathe- 
matical model for two media refraction. 

To get correct and entirely acceptable for jAiotogrammetrlc 
trangulation, the mathematical modelling must be based upon the 
actual under water position of the point rather than its refracted 
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or apparent position. The basic mathematical model for correctir 
image coordinates of under water points is illustrated in Pig. 2.1 

^d = d.h (1 - l/a)/(H - h) (2.7) 

where ^d = the correction in meters for Image point. Its 
sign is always negative (correction towards the 
photocentre), 

d = radius of Image point in the form (x + y )^ 
h = depth in meters of the underwater point at the 
time of photography, 

H = flying height in meter. 

a = ratio of tangents of angles of refrection(r ) and 
incidence (1). 

l.e. a . - 1) tan^r] "" 

where, tan r = , for vertical photograph, 

pi = index of refrectlon for ray passing from water to all 
f s= camera focal length in meters. 

If jl a 1.34 

a *[ 1.7956 + (0.7956) d^/f"^] (2.8) 

The relation between refraction corrected coordinate value 




(x' # y* ) and the photo image coordinate (x# y) is shown below 


x* = x(l - ^d/d) 

y' SB y(l + Ad/d) (2,9) 

The change in the value of equation (2,8) and its 
effect on ■<^d caused by 1® or 2® tilt is insignificant in 
aerial photograph if H is greater than 100 h. Therefore, the 
refrection ccropensated coordinates of all the vertical control 
point can be determined. To be theoretically correct, the 
unknown depths (h) of the test point should be solved as 
additional unknown in the aerotriangulation process using 
equation (2,7) through the iterative process, 

Moore (1947) developed a transparency method for water 
depth from aerial photography alone. This was one of the first 
remote sensing method. Basically the method consists of determin- 
ing the optical extinction coefficient of the water column from 
the brightness measured at identical points in two specially 
filtered (red and green) photographs, “Ihe spot depths are deter- | 
mined by means of a ratio calculator which compares the brightness 
in two band to know water colour-transmission. Results with +10% 
accuracy to a depth of about 6 m are achieved over homogeneous san 
bottoms. This technique is limited to calm sea, relatively clear 
water, bright skies, without cloud patterns and sun angle between 



2,2 AREA UNDER STUDY 


Frcsm the work of Ross^^^, Polcyn and Lyzen^^^, Jha an3 
Satya Prdcash^^^ and others it is clear that MSS band 1 reflectanc 
values could be used to predict shallow water depth* Satisfying 
above condition, coastal area of Andhra Pradesh was selected. The 
area is along the east cost of the Andhra Pradesh region between 
Madras and Kakinada. 

The area is covered by one Landsat — 4 scene. For 
reference the index map of landsat - 4 coverage is shown in Pig. 2. 
The index map is supplied by NRSA, Any scene is located by Path 
number and Row number. The area selected is located by Path 
number - Row number : 152 - 049, The area inclosed by the 
longitudes and latitudes is shown in Pig, 2,3, 

2.3 DATA COLLECTED FOR THE WORK 

The data collected for the work are as follows: 

(1) landsat data. 

(2) Topo sheets with water depth contours. 

i 

The Landsat data in the form of C 3 CT containing reflectance : 
values for band 1, 2, 3, 4# positive filmsd : 1 M) and paper pclnt| 
(1 : 1/2 M) were obtained from NRSA, Hyderabad, While selecting 





FID 2.3 AREA UNDER STUDY 



the date of pass, for which Landsat data is obtained care 
is exerci^d so that there is no cloud cover. The date of 
pass for present data is for 28th September 1983, 

The toposheets of the area were obtained from 
SURVEY OP INDIA, The area is covered by two toposheets of 
scale 1:250,000 (1 : 1/4 M). With reference to index map of 
Survey of India, the two maps are designated as 66 A & E and 
65 H, To obtain the information about availability of map of 
required scale, costs, whether restricted or not one may write 
to MAP SCALES OFFICE, HUDA CCMPIEX, TARNAKA, HYDERABAD, On 
these two toposheets itself, the two water depth contours of 
9 m and l8 m are shown. 



CHAPTER 3 


METHODOLOGY 


3.1 INTROmCTION 

To obtain the reflectance value of a ground pointy it 
is necessary to convert the longitude and latitude of that 
ground point to corresponding line number and pixel number. 

For this it is necessary to determine all the corner coordinates 
of the imagery very precisely. The corner points coordinate 
are as shown in Fig, 3.1, 

3.2 CONVERSION OF GEOGRAPHICAL COORDINATE TO CONICAL ORTHOMORPHK 
COORDINATE 

Projection of geographical coordinates of Landsat imagery, 
which is on curved surface is to be converted to a flat surface 
whose X and Y coordinates can be computed. This map projection | 
used for this purpose is the conical orthoraorphic projection 

i I 

which is used country like India to minimize projection distorti| 

The orbit of Landsat makes an angle about 9® with the 
North direction. Hence the edges of image are not parallel 
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F IG 3.1 



line no. and pixel no. 





to the longitude lines. The geographical coordinates of the 
corner points are converted to conical orthomorphlc coordinates 
with reference to chosen coordinate axis system parallel to 
longitude and latitude, 'Ihe origin Is chosen to be a point 
just outside the Imagery preferably Just below it and slightly 
towards the left so that no coordinates (X, Y) become negative 
as shown in Fig. 3.1, 

For the conversion of geographical coordinates to 
conical orthomorphlc coordinates^ latitudes and longitudes of 
the four corners of imagery should be precisely known. These 
value can be obtained from the records published by the Eros 
Data Centre, Sious Palls, South Dakota (The center supplies 
Landsat-imageries and Tapes on a world-wide basis). Pig, 3,2 
represents the four corners of the imagery. The origin O is 
chosen to be point outside the Imagery as reference point, 
let the origin O has geographical coordinates (0^ , L^) where 
0^ is latitude and is longitude. Again let (0p/ bp) be 
geographical coordinates of any point on the imagery. The 
conical orthomorphlc coordinate can be given as:- 

X = (P - m* ) sin Y 

Y » m* + X tan (•-^) 
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Fig. 3.2 Transformation of geographical 

co-ordinates to conical orthomorphic 
co-ordinates 



where 


, Y = sin (-^ ) 


Ah « 



m* 


m + 


m 


o o 


tan ta^(5 + 3 tan^ ^o) 

24 R„n| 120 


Rq = Radius of curvature of the earth (Meridianal) 

= a(l - e^) 

(1 - e^sin^0Q)^^^ 

a = major axis of the reference ellipsoid 

= 6377277,6 m (For Everest Ellipsoid adopted for India) 
b = mincar axis of the reference ellipsoid, 

* 6356075,0 tn (For Everest Ellipsoid adopted for India) 

e^ = eccentricity = 1 - h^/a.^ 

For the Everest Ellipsoid^ e => ‘ ooH " q — 


* normal to the surface at 0^ (Normal section radius) 

‘ (1 - 


P * cot 0^ 



“ W®®*’ radius of the earth (Normal section radius) 

= 

* O O 

A program to calculate the conical orthomorphic 
coordinate of a point given its longitude and latitude is given 
in the Appendix - 1. 

Here origin of the reference axis was taken just below 
the third corner point and has its cocardinate E 79® 50* 34,40*', 

N 14® 45* 13,48", The corner points 1, 2 and 3 were converted 
to conical orthomorphic coordinates. The conversion results are 
shown in Table 3,1, 

3.3 calculation OF THE CONSTANT OF THE AREA 

Transformation from conical orthomorphic coordinate to 
line number and pixel number can be made by the following 
equations. 

X « A^x + Bj^y + C^ 

Y * A 2 X + 327 + ^2 

where X « conical orthomorphic coordinate of the point s Longitude 
y =s conical orthomorphic coordinate of the point i Latitude 

X - Line number 
y * Pixe 1 number 



TABLE 3.1 


CONVERSION OP CORNER COORETNATES 


Corner 

Point 

Geographical coordinate 

Conical orthoraorphic coordinate 

Longitude 

Latitude 

Longitude 

(X) 

Latitude 

(Y) 

1. 

80.232558 

16,772421 

41552,5824008 

224769.6131003 

2. 

81.901163 

16.516816 

219743.0386210 

197273.3419803 

3. 

79.842890 

15.003363 

0,000 

27783.9174497 


Now we have six equations and six unknowns Aj,, A^. 

B 2 and C 2 . These six simultaneous equations were solved and unknowns 
were determined. 

The calculated values of the constants are:- 



s 

-17.32079300 

®1 

S 

55.04802500 

^1 

SB 

41514.855000 

^2 

SS 

-32.111587 

®2 

- 

-8.4943687 

<^2 

= 

224860.220000 


1, 


^1' ^ 2 * ^2 constants whose values can be 

determined from a minimum of three ground points. The three 
points chosen are three corners of the imagery whose (x, y) 
coordinates are already computed from early transformation 
from their geographical coordinates. The (x, y> coordinate 
(line number and pixel number) of these three points are also 
known, since they are corner of the imagery. The set of 
equations are: 

xi - hH * ®iyi + =1 

Yj - AjXi + BjYi + C2 
Xj = AjXj + BjJTj + Cj 

Yj . AjX2 + BjYj + Cj 

Xj - AjXj +8^73+ Cj 
Y3 - A2X3 + B2y3 + C3 

Here suffixes 1, 2 and 3 refer to corner points 1, 2 and 3. 

3.4 DETERMINATION OF LINE NUMBER AND PIXEL NUMBER 

Now the line number and pixel number of any point on 
the imagery can be found by solving the two simultaneous equation 


X = 

Y = A 2 X + B2y + C 2 

Here Aj^, Cji, A 2 ^ B 2 , C 2 » X and Y are kncwn and x and y 
can be calculated, A program for solving these equations 
thereby determining the line number and pixel number of any 
point on the imagery is given in Appendix — 2, 

By giving geographical coordinate in degree, minute 
and second, the program calculates the line number and pixel 
number of any point within that area. 

The scale of toposheets is 1 :250,000, The latitude and 
longitude of all the necessary points was determined precisely 
with toposheet, 

3.5 READING TFIE COMPUTER COMPATIBLE TAPE (CCT) 

The input data for the present study were taken from 
CCT. Each pixel is encoded in a byte on the CCT, Each byte 
is composed of eight binary digits (bits) which are arranged to 
represent different brightness value as binary nixnber. 

The data is stored in CCT in three files; 1st and 2nd 
files being Tape Directory and Header Record respectively. The 
tape directory identifies the contents, format of ccT etc. The 
header record givesthe content of data and discribes the format 



in which the data are recorded. Third file contains the data in 
four bands. Third file consists of 9600 records and each 
record contains 3596 bytes. Each record in third file represents 
one scan line in one band. Each scan line requires four records 
to be completely described in all four MSS bands. Ibe mode of 
data storage is Band-Inter-Leaved (BIL), In this type of data 
storage, first four records correspond to scan line number 1. 

One for each band i.e. first record gives the data for line 
number 1 in band 1, second records gives the data for line 
number 1 in band 2, third records gives the data for line number 
1 in band 3 and fourth record gives the data for line number 1 
in band 4. Fifth record to eighth record corresponds to scan 
line number 2 in all four bands, i.e, fifth record corresponds 
to band 1, sixth record corresponds to band 2 and so on. So 
in this tape, to get the data of xth line in band 4, one has 
to read '4x'th record. In each record before and after the 
data, there are zero fills. The number of zero fills varies 
from record to record. If there are n initial zeros in the 

record and if we want the reflectance value of Nth pixel, 

"til. 

then we have to read the reflectance value of (N + n) pixel 
(Byte), 


The CCT used in EEC-lO system has 9 tracks, A frame 
consists of bits 0 to 7 in which the pixel brightness is 
recorded and a ninth parity bit, A word is made of 4 bytes or 



36 bits. This is shown in Pig, 3,3 


A program for reading CCT is given in Appendix - 3, 

This program reads one scan line in one band of the image 
file. The CCT as obtained could not be read on lEC-lO system. 
This is because DEC-lO is 36 bits word machine and CCT is 
formated for 32 bits word machine. If one tries to read the 
CCT as obtained the DEC-lO system will not read the 5th, loth, 

I5th.... bytes due to lack of compatibility between the 

tape format and DEC-lO system. A dummy blank was introduced 
after 4th, 9th, l4th, .bytes increasing the record size 
from 3596 bytes to 4500 bytes. When the tape is read on DEC-lO 
system, the introduced dummy blank bytes will be skipped and 
data of all 3596 bytes will be printed. 

The required record number to be read, is found by 
multiplying the line number with 4 and adding the band number. 
All the record numbers are stored in the increasing order and 
CCT is read record by record. To get the reflectance of a 
pixel, the number of initial zeros are counted in each of the 
record and added to pixel number to get the byte number 
corresponding to that pixel and the reflectance value is read 
from CCT for that byte. 




F*9.3.3 -Cross section of q Strack CGT 


CHAPTER 4 


NUMERICAL RESULTS 


4,1 RELATION BETWEEN REFLECTANCE VALUE AND WATER DEPTH 

In water area, the light reaches back to MSS (sensor) 
as reflection from water surface, scattering directly in 
atmosphere, reflection from bottom (seabed) and reflection or 
scattering from the object inside the water. Scattering directly 
from atmosphere and surface reflection is same for same condition. 
But reflection from bottom is affected by suspended material which 
can scatter the light. The amount of light returned from the 
bottom depends upon the depth of water, its attenuation coefficient 
for the wave length and reflection coefficient of the bottom. The 
light reaching the bottom is also affected by the slope angle E' 
which can be calculated, if refractive index of sea water is known. 
Light path from sun to satellite is shown in Pig, 4,1, 

For Water depth d, the optical path length in water 

= d(l + cosec E* ) (4,1) 

If is refractive index of sea water then 

sin (90^ - E) 
sin r 


V- 


(4.2) 







where. 

E * sun elevation 

r . sln-^ ® ) 

(4.3) 


E'= 90* 

- r 

(4,4) 

Here 

E = 55* 

and |J. = 1,34 


Therefore E’ = 

64* 39* 23" 



The value of sun elevation E, is given in imagery of study 
area. When light passes through a di stance dL in water, then 
energy is attenuated to an amount dE, which is proportional to 
the energy E and the distance dL 

or dE * —K.dL.E (4,5) 


where K = attenuation coefficient and is dependent on wave 
length. The negative sign indicates loss of energy. If 
and E^ be energy of radiation for L *= 0 and L ** L respectively. 


Then, 



(4.6) 


E- = E^ , e"*^^ (4.7) 

It o 

i.e. is the energy of light when it just inters the surface. 
L is optical path length in water. Here L = AB + BC 



ca: 


L * <3(1 + cosec E* ) 


(4.8) 


Therefore E 


L 


g + cosec E* ) 


(4.9) 


Energy received by the sensor will be sum of energy 
directly coming to sensor due to atmospheric scattering and 
refracted energy coming from water E^ . E* )^ 

If E,j, be the total energy received by the sensor then 

E^ » Ed + AEq Q-Kdd + cosec E* ) (4.10) 

where E^ = energy received directly through scattering by 
sensor and Is constant for same conditions. 


A ss parameter depending upon state of atmosphere, 


or 


E~ •“ Ej = AE e 
T d o 


-Kd (1 + cosec E’ ) 


(4.11) 


If E and E be the energy (Reflectance value) received 

X y 

by sensor for water depth x and y respectively, then 


Ey - E^ 


*<3 ^ g-K (1 + cosec E* ) (x - y) 


ln( 


or 


X - y * 




-K(l + cosec E* ) 


(4.12) 


(4,13) 


If 


y * O, then 



depth X = 


(4.14) 





-K (1 + cosec E* ) 


Where E^ is energy (Reflectance value) received by MSS. 
when water depth ^ 0. 

In the above equation E^ and K are unknowns. For 
the determination of E^. E^ and K, 30 points were taken on the 
toposheets which lie on the known depth countours. Latitudes 
and longitudes of all the points were calculated precisely. 

Their line number and pixel number and finally the corresponding 
reflectance value were read from CCT which is shown in Table 4,1, 
Now we have thirty equations and three unknowns, E^. E^^ and K 
were calculated by Least Square Method, 

Solution by least square is given by the equation 

X = - (A*^ PA)*^ A*^ PL 

where A * coefficient matrix 

= matrix of partial derivative of the function with 
respect to unknowns, 

P = weight matrix. 

Here weight matrix is taken as Unit Matrix. 


L . (L^ - I^) 
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TAaLK : 4.1 

f’. <T •■.Ti. AT"’t; rAK,F-j O*' Ig'n H^A'I’ER DEPTH COHfnUR 
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TAsriE s 4.1 (....rn^'iT.'tnED) 

Are TAKEM on 9m WATER ngPTH COMTOOR 


m w* «w «w» w 

■n , 

liOS 

f.QMGTTmip'* 

wM aM m 

LTWE 

PIXEL 

REELECTAWCE 




— 

Po, 

value 

1 . 

1b, 26667 

91,53333 

4vJ4 

2655 

37 

n 

■ ■ m 

1 Is .16667 

91.33333 

576 

232? 

37 


Ib.unooo 

91.26667 

90R 

2268 

36 

■U 

15,8^667 

91,13333 

1010 

207 4 

35 


15.66667 

«l.':)0ooo 

1300 

1908 

37 

r. , 

16,66667 

«>.). 73333 

1353 

1405 

36 


16.66361 

SO, 77667 

1349 

1488 

36 


16.76600 

90.37833 

1315 

702 

35 


15,66333 

80,91556 

1321 

1750 

36 


IS. 79667 

80,67556 

1194 

124? 

37 

1 1. . 

lS,47f>78 

80,35722 

1686 

779 

37 

1 ^ 

.i * 

15,51472 

?0. 49778 

1602 

1006 

37 

13, 

IS. 36139 

90,25000 

1850 

622 

36 

lA. 
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37 



“ Matrix of computed value of x (depth) by assumed 
unknowns 

Lj^ = Observed value of x (depth) 

X = Correction Matrix for unknowns. 

The expression can be write more simpler as 
X = - (A^A)"^ A*^L 

where, X^ = Matrix of adjusted value of unknowns 
Xq * Matrix of assumed value of unknowns. 

For the assumed values of unknowns K, and correction 
matrix X was determined as explained above and correction 
matrix was added to assumed values of unknowns. This process 
is repeated till we get very less variation in the correction 
matrix. 

A program for least square method is given in Appendix 
and value of unknowns are as follows: 

» 23,61 

E^ = 62.18 

K = 0,057 m"^ 



4.2 CALCUIATION OF WATER EEPTH 


Once the unknown and K are determined, then it is 

very easy to calculate water depth at any point by equation 4,l4, 
if reflectance value of that point is given, 

4.3 DISCUSSION OF RESULTS 

There are only two depth contours on the toposheets. 

For checking the accuracy of the work, 15 points on each 
of the contour lines were taken. Latitudes and Longitudes of 
all those points were calculated. Line numbers and pixel 
numbers of those points were calculated and reflectance values 
of all those points were read from CCT, Water depths of all 
those points were calculated by equation 4,13, Ihe calculated 
depth of all thirty points are shown in Table 4,2, Then the 
difference of actual water depth taken from toposheet and 
calculated water depth is calculated for both the set of data 
9 m and 18 m water depths separately. 

Standard erroiJ (Se) and coefficient of correlation (r ) 
are calculated as: 



N 


n2 xy - (£ x) y) 

y[ nS x^ - (sx)^[ ns y^ - (S y)^]; 


r « 
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For the water depth of 9 m 

r = -0,99902 
Se s= 0,684 ra 
a - 0,706 m 

For the water depth of 18 m 

r = -0,995 

Se = 1.255 ra 

o » 1,2828 m 

The negative sign shows the reciprocal relation. Here we find 
that standard error for l8 m water depth is more than that of 
9 m water depth. It may be due to the reason that for same 
difference In reflectance values of deep and shallow water, the 
variation in water depth is more for deep water than for shallow 
water , 

From the correlation analysis, it is clear that coefficient 
of correldl:lon is about one, which shows that results obtained are 
quite satisfactory. 

4.4 IRAWING OF WATER DEPTH CONTOUR 

There are two methods for drawing contours 

(a) Direct Method 

(b) Indirect Method 




In direct method of contouring, contour to be plotted 
is actually traced on the ground. We are concerned about only 
those points which lie on the contour. In the indirect method, 
scxne suitable guide points are selected and guide points need 
not necessarily be on the contours. The guide points serve as 
a basis for the interpolation of contours. In linear inter- 
polation, it is assumed that water depth variation between two 
guide points is proportional to distance. 

In the present work, indirect method of contouring was 
done. Total area was divided into small grids by drawing lines 
parallel to the longitudes and latitudes. The interval between 
the lines drawn parallel to latitude and longitude are kept 
4 min and 2 min respectively. The latitudes and longitudes of 
all the corner points of the grids were calculated. The corres- 
ponding line numbers and pixel numbers were calculated by 
previously mentioned progran. The reflectance values of all the 
pixels corresponding to grid corners were calculated and finally 
the depth of those points were calculated. The values are 
tabulated in Table 4,3, Points of equal depths (9 m and 18 m) 
were calculated by interpolation and contours were drawn. The 
contours what we obtained, is more or less same as the contours 
of the toposheet , 
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1 o. niia 

1 a. at!enu 
1 h.7fb^7 
■> o . 2 a ^ i 
1 f>. 2^’Ouij 

itj.naaa 


’ &, iOuOu 
t b.2^Jo67 
tfa.2^3.ia 
1 o,2<nu{?C 

1 b , 1 3 3 a 3 
1 1>. 1 fiyny 
tb.iOOOO 
1b.26hb7 
16*23333 
16.20000 
16. 10667 
16.1 33 J3 
16.10000 
16.33333 
Id. 30000 
16,26667 
16.23333 
16.26000 

16.16667 
16*13333 

16.10000 
16.06667 
16.03333 
lo.yoooo 
16.33333 
16.30000 
16.26667 
16.23333 
16.20000 

16.16667 
16.133-33 

16.10000 
1 6. u6667 
16.03333 
l,6.CO0C0 
16.33333 


LPftiOiTUDE 


B1. .SOOOO 
HI .BOOOO 
HI .HOOOO 
HI .BOOOO 
61,80000 
81 , BOOOO 
81.80000 

61.73333 

81.73333 

81.73333 

81.73333 

81.73333 
81,73333 
81 .73333 
B1 ,66667 

81.66667 

81.66667 
81,66667 
8 1 ,66667 

81 .66667 

81.66667 
aJ. 60000 
81.60000 
81,60000 
81.60000 
81,60000 
81,60000 
81,60000 
81.60000 
81,60000 
81,60000 
81,60000 

81.53333 

81.53333 

81.53333 
8J.53|33 

81.53333 

81.53333 
al. 53333 
81 ,53333 

81.53333 
81,53333 
81,53333 
81.46667 


LTwE 

Mo. 


PIXEL 

r+0, 

'3127* 

3142 

3156 

3170 

3185 

3199 

3213 

3017 

3031 

3045 

3059 

3074 

3088 

3102 

2891 

2906 

2920 

2934 

2949 

2963 

2977 

2752 

2766 

2781 

2795 

2809 

2823 

2838 

2852 

2866 

2881 

2895 

2627 

2641 

2655 

2670 

2684 

2698 

2712 

2727 

2741 

2755 

2769 

2502 


reflectance 

VAL'JE 


WATfclR 

DEPTH 


20.3 
27,7 

38.3 
1,5 


14.7 

20.3 

27.7 

38.3 
2,3 
6,6 

11.8 

16.4 

20.3 

27.7 

38.3 
0.7 
3.8 

7.7 

il;f 

n:t 

18,1 

23.2 

27.7 

38.3 

2,6 

5.7 

8,2 

m 

12.7 

16.4 
20.3 

20.3 

27.7 

38.3 

2,6 
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' i , 

■ t*. 


’’ ' t-'- , 
1 7 , 

■‘I . 

''I . 

SJ, 

'■>b , 

S<*, 

'■:! . 


'I i • 

i. 

^ ’> • 

*■'' • 

fi / . 

f'lU 

7 (.1 , 

7], 

73, 

74. 
74. 
7u, 
7 7. 
7H, 
7'7, 
♦<(>, 

I * 

it 4 , 

*' 4 • 

Hh, 

W i? , 
'7'>. 
’■71, 

'■> 2 . 

43, 

44, 

45, 

4fj, 


’ i>. 
1 tt , 

V^, 
1 *> , 
1 1,1 . 
1 4. 
1 3 , 


1si\ r I T;|nf.J 


7 *> , jf’t.'Ui 

I , 27’ hft 7 
• 11,?. 3 ,i3 j 
1 n, 2 f’. 7 : 1 u 

,1 3 i 3 3 

ii'uOi.i 

• 3333 
lif'uOl; 
.4f>t-j47 

14,43333 

10.33333 

1 ». 3 ('i,Oi, 
1 u, 27657 
' 7 • 2 3 1 3 3 

U. 2«iunt; 

16.15657 

16.13333 

1 6 , A 0 6 0 U 
1 6, U7b67 
.1 6.u3 3 33 
1 O, u'tuOO 
1 6, 46p57 
,93333 
900U6 
inwO(, 
26667 
23333 
,24660 

10.16667 

16.13333 
lo.lODOU 
lb,U6p67 
Ib.u 33 !3 
1.6, 00000 
15.96667 
15.93333 
15.90600 
.15.«6667 
15.63333 
15.60606 
16.26667 
J6. 23333 
16,20000 

16.16667 

10.13333 
lo, 10600 
16.06667 
16,03333 
16.60000 
15.96667 


1 5. 

1 b ^ 
16 . 
16, 
16, 
15. 


longitude 


81.46667 
B 1 , 4 {j 6 6 7 
B 1 , 46667 
HI, 46667 
Si .4tj667 
b1 ,46667 
B1 .46667 
B 1 ,46667 

61.46667 
B 1 , 46667 
fit. 46667 
B 1 , 46667 

ai ,40000 
81 , 40000 

81.40000 

81.40000 
«1 ,40000 

81.40000 

81.40000 

81.40000 
B1 .40000 

61.40000 

81 .40000 
fal, tuOOO 

81.40000 

81.40000 

63.33333 

81.33333 
Bl. 33333 
61.33133 

81.33333 

81.33333 

81.33333 

61.33333 

81.33333 

81 .33333 

81.33333 
81.33333 

81.33333 

81.33333 

81.33333 
81.33333 
81.26667 

81.26667 

81.26667 

81.26667 

81.26667 
81.26667 

81.26667 

81.26667 

81.26667 

81.26667 


LTI’IE 
Ko , 

» •» f« » 

374 

418 

462 

506 

549 

593 

637 

6S0 

724 

76B 

R12 

855 

344 

388 

432 

475 

519 

563 

607 

650 

694 

738 

781 

825 

869 

913 

401 

445 

489 

533 

576 

620 

664 

708 

751 

795 

839 

882 

926 

970 

1014 

1057 

459 

502 

546 

590 

634 

677 

721 

765 

808 

852 


»1XEL 
14 0, 

2516 
2530 
2544 
2559 
2573 
2587 
2601 
2616 
2630 
2644 
2658 
2673 
2377 
2391 
240 5 
2419 


2462 

2476 

2490 

2505 

2519 

2533 

2547 

2562 

2266 

2280 

2294 

2308 

2322 

2337 

2351 

2365 

2379 

2394 

2408 

2422 

2436 

2450 

2465 

2479 

2155 

2169 

2183 

2197 

2211 

2226 

2240 

2254 

2268 

2282 


REFLECTAMCK 
value ■ 

43 
37 

34 

34 

31 
30 
29 
28 
28 
26 

25 
24 
61 
49 

43 
37 

35 

32 

U 
28 

27 

26 
26 

24 
61 
54 

48 

44 

37 

36 

35 

32 
29 

29 

28 
2? 

26 

25 
24 
61 
59 
56 

53 

49 
40 

|7 

36 
34 


water 

DEPTH 


5.7 

8.8 
10,9 
10,9 
13.8 

14.7 
1.6,4 
18,1 
18,1 

23.2 

27.7 

38.3 
0,3 

3.4 

5.7 

8.8 

10,2 

11.8 

12.7 

14.7 

16.4 
18,1 
20.3 
23.2 

23.2 

38.3 

2.0 

3.8 

5.3 

8.8 

9.5 

10,2 

12.7 

13.8 

16.4 
16.4 
18, 
20 , 
2'3 » 2 
27,7 
38,3 

0,3 

0,7 

1.5 

2.3 

3.4 

7,1 

8,8 

9.5 

10.9 
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^ ' '1 « 

l» i-\ f 1 ' j t% 

UnfvGlTUDE 


PIXEL 

REFLECTANCE 

WATER 

■- ' 0 



No. 

«»*» w* «w •*» »* n* 

* f 
O 1 

\rAL0£ 

DEPTH 

• ^ , 

1 ^ . *7 3 } ? j 

H1 .,26667 

896 

2297 

32 

12#7 

' ■!" m 

i b. yuuu'v' 

31.26567 

940 

2311 

30 

14.7 

■■” * 

1 *>• b^>oh 7 

1 » 26667 

983 

2325 

29 

16,4 

,1:, * * . • 

t b « 4 1 3 3 3 

81,26667 

1027 

2339 

28 

18,1 

1 ' If# 

lSi.b(‘uOi( 

81.26667 

1071 

2353 

26 

23,2 

1 / # 

1 !3 • 7 ^ D ^5 / 

81,26667 

1115 

2367 

26 

23,2 

,« S # 

1 b. M 33i 

81 .2o667 

B1 ,26667 

1158 

2382 

25 

27.7 

1 ' ’ 1 . 

1 b» /CDOU 

1202 

2396 

24 

38,3 

j ' S* 

M> . i B b fe 7 

81.20000 

603 

2072 

59 

0.7 

i ' f' • 

1b. 13133 

B1 ,20000 

647 

20b6 

55 

1,7 

1 ‘i • 

1 b. l0(>0l 

81,20000 

691 

2100 

52 

2,6 

!**''*♦ 

1 b . 0 B b 6 7 

81,20000 

735 

2114 

48 

3,8 


1 1. . t. 3 j 3 3 

81,20000 

778 

2129 

44 

5.3 

11'. 

I b.'.'Ot.or. 

81.20000 

822 

2143 

43 

5,7- 

ill* 

1 b,9fibb7 

81,20000 

866 

2157 

41 

6.6 

n 

1 b. 93333 

81.20000 

909 

2171 

38 

8*2 

1 1 3 * 

1 b.yObOt 

81.20000 

953 

2185 

36 

9,5 

j 1 

Ib.tjfibe? 

81.20000 

997 

2199 

33 

11.8 

i 1 b • 

15.83333 

81.20000 

1041 

2214 

31 

13,6 

J 1 6 , 

1 5.8f OOu 

81,20000 

1084 

2228 

30 

14.7 

n ?. 

\^,lhbfi7 

81,20000 

1128 

2242 

29 

16.4 

A 1 b . 

lb. 73333 

81.20000 

1172 

2256 

28 

18,1 

11 y. 

1 b.70uOO 

81,20000 

1216 

2270 

27 

20,3 

i?i*. 

15,66667 

81,20000 

1259 

2284 

26 

23.2 

i?i. 

15.63333 

81,20000 

1303 

2299 

24 

38,3 

172, 

15.96667 

81.13333 

879 

2032 

59 

0.7 

I7.i, 

15.93333 

81,13333 

923 

2046 

51 

2,9 

174. 

15,90006 

81.13333 

967 

2060 

48 

3.8 

1?5, 

15,86667 

81,13333 

1010 

2074 

43 

5,7 


15.83333 

81,13333 

1054 

2088 

37 

8*8 

177. 

15.80000 

81,13333 

1098 

2102 

36 

9,5 

1?H. 

15.76667 

81.13333 

1142 

2116 

u 

12,7 

17 y , 

15,73333 

81,13333 

1185 

2131 

29 


117. 

15,70000 

81,13333 

1229 

2145 

28 

18,1 

1 1 1 , 

15.66567 

81.13333 

1273 

2159 

27 

20.3 

1 u. 

15.63333 

81.13333 

1316 

2173 

25 

27.7 

11 j. 

15,60000 

81,13333 

1360 

2187 

24 

38,3 

114. 

1 5.90000 

81.06667 

980 

1934 

58 


1 ly. 

15.86667 

81,06667 

1024 

1949 

52 

2 ,6 

1 IB. 

15.83333 

81,06667 

1068 

1963 

46 

4,5 

13 7. 

15.80GQU 

81.06667 

1111 

1977 

42 

6,2 

i i y . 

15,76667 

81.06667 

1155 

1991 

39 

7,7 

1 17, 

1 5. 73333 

81.06667 

1199 

2005 

36 

9,5 

1 4 1 ) . 

15,70000 

81.06667 

1243 

2019 

32 


141, 

147. 

1 4 i , 

15.66667 

81,06667 

1286 

2033 

30 

14,7 

15.63333 

81,06667 

1330 

2047 

28 

18,1 

15.60000 

81.06667 

1374 

2062 

26 

23,2 

14 4. 

1 5,56667 

81.06667 

14^7 

2076 

25 

27,7 

i 1 5 • 
i4b* 

1 47 

^ 15*53333 

81,06667 

1461 

2090 

24 

38,3 

15.73333 

81,00000 

1212 

1880 

46 

4,5 

15.70000 
15*6666 J 

81,00000 

1256 

1894 

40 

7,1 

4 " f • 

4 4 M . 

81.00000 

1300 

1908 

36 

9,5 
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T' , 
tM, 
.1 ^ • 
!H .. 
IH' , 
I' I. 

I ' ' # ♦ 

1*’ }. 

l^'C » 

1^7. 

1 (t H , 

H.'i, 

i » 

1 M . 

i 

\u» 

I’f'K 
1 'fs. 
llh, 

r? ;. 

* '7H, 

I 

M'-i. 
1 H J . 

i'u. 
14 4, 
i H b • 
r!n. 

i. 'Ui • 
1 '^'4. 
V’*u, 

rn. 

i . 

j '•* 1 • 

V'’4. 

I 

1 6 , 

,i 7 , 

i-7y. 


I *> . 7 

I 'j » 7 £■• i) 0 y 
" 7 

’ 7 . 6 4 i 4 3 

^ • f'* B y i ) iJ 

1 !) » h b 7 

1 "j . b ’M 3 3 
1 :5.7« Ot-O 
1 b,()fco67 

1 j • ()O0fiy 
1 '> ♦ 7 

1 *3 . S 3 j 3 3 
is.brt uou 
lb. 7CU0U 
1 b.t>f>b67 
1 b. 63333 
tb.tnoot) 
lb. 56667 
1b. 53333 
lb,5fJCCC 
15,66667 
15.83333 
15,80006 
15.76667 
15,73 3.33 
15.70000 
15,66667 
15.63333 
15.60000 
15.56667 
15.53333 
15,50000 

15.46667 
15.43333 
,15.40000 
15.86667 
15.83333 
15,80000 
15.76667 
15.73333 
15.70000 
15,66667 
1 5.63333 
15.60000 
15.56o67 
15.53333 
15,50000 

15.46667 
15.43333 


*^’*****''***fw<«»ii*««** 

6''*‘ViC?,lTUDE 


i3 1,00000 
81 .00000 
81.00000 

80.93333 

50.93333 

80.93333 
«'N93333 
BO. 9333 3 

80.93333 

80.86667 

80.86667 

80.86667 
80,86667 
80,86667 
80,86667 
80,86667 
80.80000 
80,80000 
80.80000 
80.80000 
80.80000 
80,80000 
80,80000 
80,73333 

80.73333 

80.73333 
80.73333 
ao.p333 

80.73333 

80.73333 

80.73333 

80.73333 

80.73333 
80.73331 

80.73333 

80.73333 

80.73333 

80.73333 
80,66667 

80.66667 

80.66667 

80.66667 
80,66667 
80,66667 
80,66667 

80.66667 

80.66667 

80.66667 

80.66667 

80.66667 
80,6666? 
80.66667 


i«it^ 


LINE PIXEL Reflectance 
__No. No, value 

1^22 28 

13o7 1936 25 

1431 1950 24 

1269 1768 44 

III! UP 39 

1357 1796 31 

1401 1610 26 

1444 1824 25 

1488 1839 24 

1283 1643 48 

1327 1657 46 

1370 1671 43 

1414 1685 37 

1458 1699 27 

1501 1713 2S 

1545 1727 24 

1296 1517 48 

1340 1531 44 

1384 1545 38 

1427 1559 32 

1471 1573 27 

1515 1587 25 

1559 1601 24 

1091 1321 58 

1135 1335 52 

1178 1349 48 

1222 1363 43 

1266 1377 41 

1310 1391 39 

1353 1405 37 

1397 1419 36 

1441 1433 32 

1484 1447 .30 

1528 1462 28 

1572 1476 27 

1616 1490 26 

1659 1504 25 

1703 1518 24 

1104 1196 58 

1148 1210 51 

1192 1224 48 

1235 1238 44 

1279 1252 39 

1323 1266 38 

1367 1280 36 

1410 1294 35 

1454 1308 34 

1498 1322 31 

1542 1336 29 

1585 1350 28 

1629 1364 27 

1673 1378 26 


water 

DEPTH 

27.7 
38.3 

5,3 

7.7 

13.8 

23.2 
27.7 

38.3 

3.8 


8,8 

20.3 
27.7 

38.3 


8,2 

12,7 

20.1 

U:l 

3,8 

5.7 

§'•§ 

7.7 

8.8 

9.5 

12.7 

14.7 
18. 1 
20,3 

23.2 

27.7 

38.3 
1,0 
2.9 
3.8 
5.3 
7.7 
8.2 

9.5 
10.2 
10,9 
13,6 

16.4 
18,1 

iM 

23,2 
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‘-ri. 
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> 'W« («« jw» M Ml 
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i " i , 

1 ! • ijk , 

2 ' ; ij . 

4 ”' t' • 

4 

<ii a. 

21 1 . 
212. 
21 i, 

214, 

215, 

216, 
21 ?. 

ils; 

ill* 

ill' 

♦ 5 * 

227* 

228. 

ill: 
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215. 

23b. 

237. • 
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25'J. 
210, 

2 - 1 1 . 

2-U, 

2 H. 

2 4 ‘1 » 

2 1 j • 

2 1 6 . 
217. 
2l.i. 
2 13. 
2 5 ■.} . 
2S i. 


l.Hril’urij;: 


15.^0jD0 
15 . 3655 / 
15.33333 
1 5. 3556 ? 

1 5. ^333i 
1 b. 3 POOJ 
is. 75 o57 
15 . 7 3 5 3 3 
1 5,?r>,>nj 
1 5* 65b5? 

15.53333 

15.50003 
1 5.'5§b57 
1 5,53333' 
15.50003 
15,46067 

15.43333 
^ 40003 


1 5 


n-.m 


67 

33 




i^OHGl'fUyZ] bil'MS'- ' 

mlm'm 'm 'm ’ « ]#* |«. ]*» \m j** ^ 


1 

15. -71333 
15.-70005 
15.'6|667 
15,‘5|3|| 

11 : 15 ^ 1 ? 


15.40000 
1 5.-3656'7 

mini 

1 5 . -26667 
15.23333 

is.-aoooo 

15. '7 6567 

15.73333 
15. '7 0 00 0 
15. 56567 
15.53313 
1 5 .- a 0.>05 
15.56557 
15, 53333 
16»'6230D 
t 5 , 46 b6 7 
1 5 • 1 3 3 3 3 


30.65657 

30.65567 

30,6556? 

30.63000 
30. 63033’ 

30.63033 
aO. 63003 
32.60003 
aO.'Soooo 

80.63003 

80.63003 
82.f22i>5 

80.63033 
80.63000' 
80,63003' 

80.63003 

80.63000 
50,63300 
sO. 63000 

:Um 

i' 

80. § 333-3 
80,63333 


I7l5 


tkilirt 

80.63333 
8 0.6333 3 
80.63333 
80,45567 
30,45667 
80.46557 
80,456d7 

53.45657 

^o:A^i§7 

80:45667 
80. ,45567 


1834 

Ml? 

11 51 

IW 

iW. 

l?9?' 

1336' 

13'3 0 


145 7 


159^ 
16 4 2- 
15B5' 

fill 

1174 
I - - 

I'l^OS' 
1349 
139 3- 
14 3-7 


nil 


1599 
13 43- 
178 7 


nil 

iiii 

nil 


41 1 
1494 
163 a 

nil 


Piicei, 

»'**^ «•>'** iW -Mir' 

13 92 

1 40b 

1420 
1 070 
1084 
1098 
1112 
1125 

1 If a 

1154 
1 16 b 
1^2 
1 19& 

1 210 
1224 

1238 

nil 

’i§3 




lOOl 

lots 

1029 

1043 

1084 

1098 

Ilk 

iisS 


1 |H Ull 


1182 

U9b 

1210 

847 

861 

875 

889 

9 O 3 

917 

93l 


lOOl 


***'«»« |»i* '»•]**:***• *«»W‘-s 
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i/R>tOp’ " 

• j«» m'm wW «,]** ««> M* *•■}«*', 

25' 

2 5^ 

24 

59 

48 

43 - 

4 3 
37 

ll' 

3 4 
3 3- 

II 

k 

27 

25' 

2 5- 
2'5; 

24 

51. 

4i 
4! 

39' 

?l 

33. 

7'3. 

2 I 

23 
27 
27 
2'5' 

23 : 

2'5. 

55' 

54 
43 
f5' 

43- 

Vi 

3 
3 

11 : 

30 


' •» kM * 

. P.a:PrH 

2 '8 ..2 
I'f 

27.7 

1.7 
3-. 8 : 

5.7 
7,'t 


1 0 .-9 

1 t ,8 
1 2!, 7 

1 4 .7 
18.1 
io.-i' 

2 0.3 

2.3., -2 

2 3w- 2 
2 7,7 
33,-3^ 
l‘,'2' 
3 ',3 

;5;1 

1 1*.8 
1'2!.7 
1 3r.8 
1 4.7 
I 5'.4 
1 S', 4 

1 8 .'t 
2,>,3 
20.3' 

2 1':! 

2 7.7 

3 8.3 
L •. S 

2:,0 

3'*'i 

4,5 
5 . .7 
7,1 
3 .8 
t'0,'2 
I ?;*7 
1 It . 3 

1 3t.'8 
1 4.7 
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15 
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15 
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No. 
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1638 
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26 

25 

ll 


water 

DEPTH 


48 

48 

46 

44 

42 

i? 
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26 
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32 
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Fig. 4,4: Map of Coastal Area 

* -#^tand Water having suspended materials 
^ .AClear Water 



CHAPTER 5 


CONCISIONS, LIMITATIONS AND 
RECOMMENDATIONS FOR FUTURE WORK 


5.1 CONCLUSIONS 

The following conclusions are made after studying the 
Landsat data of the study area. 

Bathymetry and coastal georaorphology can now be obtained 
from satellite data using Remote Sensing techniques specially In 
Band 1 (Landsat 4 and 5) where water is shallow and clear. The 
information obtained by this Remote Sensing technique Is much 
time saving and cost effective compared to other conventional 
method like sounding machine, echo sounding, sonar, visible 
spectral range aerial photograj^ic techniques. The scope of 
this technique becomes more, particularly in those area where no 
recent hydrograjAilc surveys are available. This project Introduce 
the newly developed technique of ‘Bathymetry* i,e, delineation of I 
coastal feature, water depth, water having suspended materials, eti 
with the help of spectral reflectance by sensor, | 

' i 
1 

The scope of present project enables us to obtain the j 

pareplanned position and to improve the internal relative position 
of the coastal features, in detecting features missed by earlier I 



surveyoJ^®' in imprcwlng <^lineatlon of features and to demarcate 
high water line and low water line. 

Remote Sensing technique of Bathymetry is highly useful 
in areas of mobile seabed which are changing frequently like 
sandspits. at mouth of rivers^ newly formed banks. Islands to 
determined the rate of change of features and to update the 
information. It further helps in the study ef coastal erosion^ 
bank , channel, ground , submerged feature etc. The delineation 
of features by remote sensing is not only time saving and co^ 
effective but also helps in developing new ports and sea routes. 

in conclusion, it may be stated that present work helps 
in the flel^ of hydrography as follows: 

To locate the high water line and low water line. 

To get water depth with good accuracy. 

To detect the character of coast and foreshore. 

To detect uncharted, mispositioned and submerged features 
which are hazardous to navigators. 

To define base line for territorial boundaries. 

To locate newly formed islands, banks, shoals etc. 

To show frequent change in ports, harbours, coastal areae' 
of which Itroedlate survey is not possible. 


(a) 

(b) 

(c) 

(d) 

(e) 

(f ) 

(g) 
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5.2 LIMITATIONS 

The followings are requirements of this study: 

1, The area should be cloud free, 

2, Water should not have suspended materials because it can 
alter the reflectance value, 

3, It should be shallow water, 

4, There should not be change in refractive index of sea 

I 

water. For our work, it is assumed that refractive index I 
of sea water is same all over the sea, 

5, The alignment of canals, drains, streams can be well 
interpreted but depth of water in canal, streams can not 
be found from Landsat Data, 

6, Landsat data is found useful for regional study of large 
water bodies rather than small areas, 

7, In general, ponds larger than 1 hectare and streams wider 
than 20 m are seen easily by Landsat imageries but 
lndlvl<Aial objects less than 70 m can not be distinguished ; 
easily by Landsat data, 

5.3 RECOMMENDATION FOR FPTURE WORK i 

1, The work can be further extended so that water depth can 
be predicted for deeper water, 

2, If water depth is calculated keeping in mind that refrectiv€ 



index of sea water is changing from place to place, we 
can get better result. 

In this project, contouring was done for clear water. 

One may extend this for water having suspended materials. 
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APPENDIX I 1 


*‘*1LE NAMF ; man, FOR 
CnMNfERT 

Tci Its conical DSTHOML„- .._^ . 

the COf^S'^.VNTS are VAUt5 FOR SfUDY ArFA OMLY, 


"nn^RAH Tn conn^ert the geographical COORPIMAXE 
' PHJg gpROINA_ X S, Y 


\r 


IMPbTCXT DOURLE PRECISlON( A-H,K-E) 

nj'^ENSTOM LAT C iSO) ,lONGC 150) , LATRADC iBO) 

I.A«Gna79,842R899 . 

l.AT0si4 7,53744l22»3, 1415921/180,0 

As5i77277,6 

ELLIP=l,0/300,8 

eli.ip»§£lip#ei,lip ^ 

SlNsDSlNCliATO) 

SINSOssSINfSiN 

DENOMaC1.0-(ELt.lP»SIflS®)) 

NOaA/CDSORTCDENOMn 

RO»NO*CCl,0-El,tilP5/DENOM) 

P=NO*pCOS(LATO)/6siN(LAfO)) 

RMsDBQRTCROfNO) 

READC21,»)II 

DO 100 lati.II • 

GAMMAaGAMMA43,14lS92t/iif»0 "'■> ■ 

M=RM*C(LAXRAD?I)-1,AT0)3 

MlsM*^l3/C6.0’t'RO*NO) 

TANaDSlHaATOJ/DCOSaAfO) 

M2=«l*Cf M»TAN3/(4,0»fiO)) 

M11sM4-M{ + M2 
X«CP-MU)YDSIHCGAMMA) 

GAMMAbGAMMA/2,0 

ysMlltX»CDSINCGAM«A)/DC0S(GA«MA3) 

»^RIXE(23,10) LATtI),LONGCI), X,Y 

CO^ITtwt^F 

F0RMAIC5X,ri0,6,2X,F10,6.2X,'X*r,F23a2,5X,'Y»',F23,15/) 

STOP 

END 
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APPENDIX ! 2 


i^!iRl"/iir^?flI&i?j|^SfiTlr^°"^*“^^ |i?ir 

COORDINATE, IM^GKRY AFTER CONVERSION TO CONiCAIi OrIhDMOR 

Vyir' nnMc-^mw€<tm^ m -« .: 


T a* rk » m> M’W ^ .. 


DlSFNSfL*^?yS^fcS?®CISXOM(A-H,K-E) 

|£p?Jk£il?g2glR 

SIN 8 qS|iH»SIN^ 

^gipMat , 0 **i£i:^jp;<esjfisQ 

|Sjg||g??i5WS?)''^«''«''*T0)) 

2§ 

"AT?|p6?|§fg}(55Jg^g5gifBW|ASI»!!<SICa 

siMiA»iiOiCDr^sii ■ * 

"fRJjMlATgROtljJiATO) ”•*" 

^lS^Jfl/ 5 ®tO*RO»N 03 

TAMaOSlNCLArO)/DCOSCLAT05 

I? 

8 i» 54 , 56466 | 733 $l 859 

Ai=^ 7 l®llllfllflHp 57980 

9250781 7511000 


''"'Tp I x«|ff |*§ ^ ^ ^ ^ 2 ^ 


END ■ ' # «.!« 
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APPENDIX I 3 


C 

r 

? 

r 

f 


”””i*-!rMjrt’c??7DR*”***"**’*’'***********"*»*«*»»»”»** 

THTs'’p»g|8j8'^SgA8g iSS^IrWBn nr rr% 

^ rnM\ff?:R.s|H5 of CCT FOpL? rSin vlSlLJI VAMD AFTEP 

opiStu'liiTSsiJolviiFlJfisPi”’'”'''"!'* 

np^NCUJV’ItjcSojD^ncIs'DSK'' 

2PSNfu?!iI=§5.&syxcl=:pir 

"C C«£^ 


li 


60 


71 

101 

1.0 


READC20)OCT 
DO 1200 Isl,U25 
NOaNO+4 

BT1»0CTCI)/2**28 
f|2aC0CTf |1-BTJ»2»*28)/2#»20 


-RECORD SliEwll25,nE: 


"oRMsIcIfl 

CONflNUE 
S.AIfL CPHE 


2«2jt20)/2*tefl2 
‘‘f3*2l»12)/2*f4 


w ri ■#«# w 'V, 111 C-, 

Ifis! , ' ■■iSC'^v ' 


';OOfl*'DOKP',RECORO »l*E.uas,DF 


END 

®ySgOU|lNE cone 

c^SSgSg^gJx"”'"" 

RIwino S3 
2S 19i,I=lf225 
gEADC53, 10) CNATCK 
Si^S^IL 

gEAD(S3,lO) (WATfK 
readcsb'ioKmatck 

loUUmnr^ 

CONTINUE 

igPHAXCSIS) 

RETURN 

Em 


,Ka 

:Kili 

-Ka2i;2 


' ■«?’ I 
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APPENDIX I 4 

45 tntt t ’( 54 : *♦*>}: Kf 4 e ft #♦!>!♦** 

FjT,t hA^-IF t r.EAST,FOR 

f^KnGRAM tq pT^n THE VALUE OF CONSTANT EO, E<f ^ Vci m LEAST 
riTTCH IS USED IN DEPTH EORMULA 

(UMEUSryW VXf l<))fX(10),AT(^,lO) ,ACtO,3)r ALClD),Rt CtO),R2flO) 

niMKMSlpV XI CIO) : 

npt.'HCU'''|'»’s 22 ,DEVlCEa' 0 SK%riLE='FOR 22 ,OAT',ACCEsS*'SEQnU') 

UPEHC 0 MiT= 7 i, 0 EVTCEa'DSK',FIT,E 3 'F 0 P 2 f.DAT^lACCEls»''SKQTN') 

X 5 - 2.267 

VOSS'? 

V 0 al 7 

CKaOjOVS 

DO m Ksl,12 

VpaVD ' 

CfCaCK 
VO»¥0 

DO 10 1*1,10 


« * * ^ n w W \ f \ •» / - ' w ^ r \ 'W w "W 

T^PE»*VX j I) *X S I L Rt C I ) f li 2 Cf;|:t'‘^-., V 
continue 
dqSo a®l,10 

ATCl,a)»Ri CJ) 

I w ? 

DO '40 Jsl.io 
ATCl, J)*R2CJ) 

D0^45 0*1,10 
ATCI, J)*R3CJ> 

Ja J 

DO 31 1*1, 10 
ACI,J)«Rlh) 

Js 2 

?■ 

0»3 

Dtt.,:4i: la. 


DO' 41 I« 
A'ri,\ 


0)*R 




0»:3 

.R! 


m 


m 


; ■' „ 


... 

La 3 
MslO 

M'ig 3 ..... ' ' ■ 

CALL MULTCAT,A,B,L,M,N) 
M a 3 

^ALL INVFRTCR,N) 


A(i? 5 ^)fn!ioj: 33 }: 3 ) 


.Rjei), 

w’si » |o 5 


1 * 1 , 10 ) 
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'■’5:1 0 

CAT.i, ;^ltLTCB,AT,C,lLifM,N5 

fj S 1 

■’stu 

'.is I 

r.tKt.u Mlft,TCC,At,,DrL,M,N) 
Oia-nt 

02-«n2 

Hia-n;? 

CRarK+ni 
Vj)a\/nH(.D2 
V(}aVOfD3 



flfPE*,Oi,D2,nS 
CONTINUE 
STOP 
END 

subroutine for MULTlPLlCAftON OF MATRIX 
SUBROUTINE MULTCR,S,T-L,H-fl) 

DIMENSION HCL,MJ,§Cfi,ft);Tll,,N) 

DO 5S I»1,L 
DO 55 Jssl> 


i' 


SUMsO.O 
DO 46 Kst.N 

8UMJsUMtRh,K)*SCK,a) 
TCl,J)aSUM 

Rl 

END 


J*' * 



■.]U' 


DO 113 _ 
lND|X(I)i 
AMAX=-1,0 

,115,117 


TCOL*! 

AMAXwTEMP 
COMTINUE 

IFCAMAX)118,118,119 
INDEX (tCOL) at 
PIVOTssAflCOOflCOrO 

AhcOU,lCOL)*!iO 

DIV0T»1,0/PIV0T 
DO 110 Jal,N 

AfICOI*f.|)aAClCOL,a)»PlVOT 

If cl-j 

TEHPPJ 

Dg ill J*|>M 

clflftloE 

00 TO 114 
return 

END 


»i # i 

I»n 22 , 12 i ,122 
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C 


APPENDIX f S 


C 

n 

c 

r 

e 


1200 


C 

C 


ILK NAf»!F: ; part, FOR I 

BUT THE i 

nurpur wtlIi be only a part of it. 

l'p?«?SSlTi5i;jilkgU6iP5''^’-BI4,BYTEC4000) 

OPE(vlCUNlTo50,OEVlCEa'l)SK',FILE«'PART'3 

Sl2fi*U25, 

READC20)0CT 

JsJ 

DO 1200 Tat. 1125 
8TlsOCTCri/2*«28 

RYTECOJaRTl; jaj+i 

ilTECJ)l8|4faJj+l*^®*®”*2**20.BT3*24»12)/2*»4 

coiriiuE r V ^ * 

SI$E»lti25, 

TO WRITE THE. DATA 

WR I TE C 5 0 , ♦ ) C BYTE CK ) , fUm) 

END ^ 
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APPENDIX 


file name s density, for f 

THIS IS AN interactive PROGRAM FOR MAPPING THE 

based on the density Slicing method I 

PKOGRAM TS FOR effective FOR 10 CLASSES 

there rs HO LIMIT ON THE LENGTH OF THE INPUT SWT ; 

NUMBER^DF PIXELS PER LINE SHOULD NOf EXCEED 960 




INTEGER I 
INTEGER C 
openconit 


¥2?i i 


-LONLIMCIO) 


'MHPl 


!!D§Kt f FliiE^XllAPi 


un. T ji.vj'j® 

OPENCUNiTalO, DEVICE! 
OpENCUNiTsli, DEVICE! 
OPENCUNITal 2 |.DEVICE! 


^ UH4. t -413 w yoR , r 4.M6» AHry T 

OPEN(UNITa6,DEVlCEa'DSK',|'ILE»TMi^Pl»} 
OPEN C UNITS,?; DgVlCE^'DSK'JriSE «%P2^ 
OPENCUNiT»|;D|viCE='DSK';FlSi«^MApP) 
OPEN CUNITaP; DEVICE® 'DSK^-'nlilss^MApPS 
OPEN? UNlTal 6, DEVICE® 'DSK*, FILE* 

OPEN CUNlT=li,DEVICE*^DSK''JptLi« 'MAPS' 
OPEN CUNITal2rDEVICEa'DSK'; FILE® 'MAP?' 
OPEN CUNITal 3, DEVICB*'DSK'J FILE* 'MAPS'! 

INTERACTION BEGINS 


Ty PE i 0 

FORMATC' ALL TERMINAL INPUT If FORMAT FRiE'iil'/. 

NUMBER FOR LkfRs!; 

TfPE 20 

fOSS^El^ArnSLi!? CHARACTERS FOR REPRESENTATION 

accept |SMCGREirCl 5 ,I»i,NCLASS) 

FORMAT CiOAl) 

level and interval of each class ' 

FROM TERMINAL < 

DO , 30 I»I,NCLASS 
TYPE 40,1 

Jg?g5L'li'§BiS??LiLISIv£?Sr"'* i»i«!->v*b',x2, 

’LONLIil|l»CONTORCl)-rNTVLCI) 

TYPE 56 


all terminal^inpot 


FORMA 






tCHMATCSX.'TYPE IN THE NUMBER OF OUPUT FILES') 
ACCEPT *.NOUTFL 


INTERACTION ENDS 


?s 



TO IWITIAOISE the IHTERmIdiXte^’maIrI jc^f I 

DO 65 Ial,NPIX 
INTMATCDa' * 

TO READ AHD CLASSIFY THE 

iLoopsso 

DO 70 lal^NFIX 
D0".70^a»i;ilCLASS„ 

CONTINUE 

intmat to output iii## 
for getttfi# tD« line printer »ap, 

MU^?lT=f» ' 

Ms?l . ■ ■ ■ 

NsM + il9 




P.LT. INLINE) SOTO 60 


